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Forward 


The high-lift calculations for this paper were made on the NASA Origin 
Super Computer in the year 2000 to 2004 time frame. Other obligations 
prevented publishing the results immediately. At the time, the NASA Origin 
Computer was the most powerful computer at NASA. Since then, other 
researchers have made high-lift calculations on far more powerful computers, 
and have advanced the art and science of numerically predicting high lift flows. 
Still, the high-lift computations of this paper compare to wind-tunnel data at 
least as well as many of the calculations published at a later date. In addition, 
many of this study's findings have been corroborated by later research (examples 
of which will be cited in the discussion). Hence, these calculations have stood the 
test of time. 

This paper is published now to document the work, since it contains a 
broader analysis (particularly of the flow field) than many other papers, and 
since it has stood the test of time. Thus this paper adds to the volume of work 
that sheds light onto, and gives confidence in, the basic technology for solving 
these types of high-lift flows. 
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Abstract 


Solutions were obtained with the Navier-Stokes CFD code TLNS3D to predict the 
flow about The NASA Trapezoidal Wing, a high-lift wing composed of three elements: 
the main-wing element, a deployed leading-edge slat, and a deployed tr ailing-edge flap. 
Turbulence was modeled by the Spalart-Allmaras one-equation turbulence model. One 
case with massive separation was repeated using Menter's two-equation SST k-co 
turbulence model in an attempt to improve the agreement with experiment. The 
investigation was conducted at a free stream Mach number of 0.2, and at angles of attack 
ranging from 10.004° to 34.858°. The Reynolds number based on the mean aerodynamic 
chord of the wing was 4.3xl0 6 . Compared to experiment, the numerical procedure 
predicted the surface pressures very well at angles of attack in the linear range of the lift. 
However, computed maximum lift was 5% low. Drag was mainly under predicted. The 
procedure correctly predicted several well-known trends and features of high-lift flows, 
such as off-body separation. The two turbulence models yielded significantly different 
solutions for the repeated case. 



Introduction 


As computational fluid dynamics (CFD) has matured, it has become 
increasingly useful in aircraft design (refs. 1 to 3). For instance, CFD is now 
widely employed in designing two-dimensional airfoils for cruise. Significant 
progress has also been made in the challenging area of drag prediction of aircraft 
configurations (refs. 4 to 6). Progress is even being made in modeling the 
complex aerodynamics of fully three-dimensional wings configured for the high 
lift that is demanded at takeoff and landing. 

When configured for high lift, wings typical of those on current transports 
usually have multiple elements and are at relatively high angles of attack. As a 
result, strong, mutually interacting circulation regions exist around each of the 
elements. Interaction between the boundary layers and wakes of the elements, 
the possibility of transonic flow with shocks on the leading-edge slat, boundary 
layer transition on each element, and the influence of three-dimensional effects, 
all complicate the flow. The possibilities of laminar separation bubbles on the 
forward elements, massive separation at the higher angles of attack, and 
relaminarization of the flow in regions of strong acceleration add even more 
complexity. 

Because of the difficulties of modeling these intricacies, early numerical 
studies of such high-lift problems were limited to two-dimensional flows. (Three- 
dimensional flows were usually handled by semi-empirical techniques. See refs. 

7 to 9.) The two-dimensional numerical calculations often over predicted the 
max i mum lift compared to wind-tunnel measurements. As advances in grid 
generation, numerical solution techniques, and computer capacity enabled the 
analysis of high-lift flows in three dimensions, three-dimensional Navier-Stokes 
investigations were conducted initially for un-swept wings mounted on wind- 
tunnel walls. (See refs. 10 to 12.) These studies indicated that the over prediction 
of lift indicated by the two-dimensional solutions was at least partly due to an 
experimental loss of lift due to three-dimensional flow, which developed at the 
juncture of the wing and the wind-tunnel walls during the two-dimensional 
tests. 
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As a result, a series of joint computational and experimental investigations 
of a high-lift wing was initiated for code validation and development. The 
configuration chosen for the project had a swept trapezoidal wing mounted on a 
body pod. It has become known as the NASA Trapezoidal Wing. The wing was 
composed of three elements: the main-wing element, a deployed leading-edge 
slat, and a deployed trailing-edge flap. The investigations were conducted at a 
free stream Mach number of 0.2. Experimental measurements included force and 
moment data, surface pressure data, off-body data, boundary-layer transition 
data, pressure sensitive paint data, PIV laser data, and acoustic data. These data 
were taken during several entries in the NASA Langley 14x22 Wind Tunnel, and 
one entry in the NASA Ames 12-Foot Wind Tunnel. See references 13 to 18. 
Reference 18 gives an account of the various investigations and the primary data 
sought for each one. 

Refs. 13 and 16 report the comparison of numerical calculations with 
aerodynamic force-and-moment and surface pressure data from mainly the 
NASA Ames 12-Foot Wind Tunnel test. These calculations were made with the 
implicit finite-volume Navier-Stokes code OVERFLOW, and with turbulent 
viscosity simulated with the Spalart-Allmaras turbulence model. Calculations 
were made for both free-air and in-tunnel boundary conditions. There was good 
agreement between the calculated and experimental lift coefficient values at low 
angles of attack, but large discrepancies between them were observed near the 
angle of attack for maximum lift. These results also showed that for the most 
accurate simulation of the wind-tunnel data with large tunnel-wall interference, 
the numerical calculations needed to simulate the wind-tunnel walls. Reference 
19 contains a comprehensive listing and review of CFD subsonic high lift work 
between the years 1988 and 2002. 

More recently (circa 2010), the American Institute of Aeronautics and 
Astronautics (AIAA) held an international workshop, AIAA HiLiftPW-1, on 
calculating the flow for the NASA Trapezoidal Wing. Compared to the previous 
calculations, many of those presented at the workshop showed significant 
improvements on the agreement with the wind tunnel data. They illustrated the 
importance of strategically refining the grid, modeling the attachment brackets 
for the various elements of the wing, modeling boundary layer transition, and 
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including the cross-derivative diffusion terms for high-lift Navier-Stokes 
calculations. See references 18, and 20 to 22 for a summary of the results. 

The present paper presents the comparison of numerical calculations 1 * * * * with 
the original set of the NASA Langley 14x22 Wind Tunnel data, data that was 
corrected for wall-interference effects using classical methods and is considered 
largely wall-interference free (refs. 13 and 23). The numerical computations for 
the present paper were made with the finite-volume Navier-Stokes code 
TLNS3D (refs. 24 to 28). Turbulence was modeled by the Sparlart-Allmaras one- 
equation turbulence model (ref. 29). For comparison, one solution at a high angle 
of attack was repeated with Menter's two-equation k-oo Shear-Stress-Transport 
(SST) turbulence model (ref. 30) since it has predicted more accurate results in 
separated regions for several configurations (refs. 28 and 30). The numerical 
calculations for the present paper were made at a free stream Mach number of 
0.2, and at eight angles of attack ranging from 10.004° to 34.858°. The Reynolds 
number based on the mean aerodynamic chord of the wing was 4.3xl0 6 . 

Comparisons are made between the numerical and experimental results for 
force and moment data, and for the surface pressure coefficients. The results 
obtained with the two turbulence models are also compared for one case. In 
addition, trends of the numerical results are compared with other investigations. 
Since the purpose of the paper was to evaluate the "off-the-shelf" performance of 
the code and turbulence models, no attempt was made to calibrate the code or 
turbulence models by adjusting parameters to improve the agreements between 


1 These calculations were made in the 2000 to 2004 time frame. Obligations 

on more pressing research prevented publishing them at that time. This paper is 

published now to document the work, since it contains a broader analysis 

(particularly of the flow field) than many other papers, and since many of its 
findings have been corroborated by later research (examples of which will be 

cited in the discussion). Thus this paper adds to the volume of work that sheds 
light onto, and gives confidence in, the basic technology for solving these types 

of high-lift flows. 
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the numerical predictions and the experimental data. Finally, the general 
characteristics of the flow field are discussed. 


Nomenclature 


Unless otherwise noted, all variables are normalized by appropriate 
combinations of the free-stream parameters, the reference length, c , and the 
reference area, S. 

b / 2 wing semi-span, inches 

L 

C, lift coefficient, 

qooS 


C D 


c 


m 


drag coefficient. 


D 

qooS 


pitching moment coefficient. 


m 


qooS c 


C P 

CFD 

constrained streamlines 


c 

Qocal 


D 

8s 

gf 

h s 

IGES 

L 

M 

m 


n — n 

pressure coefficient, — — — 

computational fluid dynamics 
computed streamlines that have been 
numerically constrained to a plane 
mean aerodynamic chord, inches 
local chord of the wing in the cruise 
(stowed) configuration, inches 
drag, pounds 

slat gap, inches (see table 1) 
flap gap, inches (see table 1) 
slat height, inches (see table 1) 

Initial Graphics Exchange Specification 

format 

lift, pounds 

Mach number 

pitching moment, inch pounds 
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Of 

p 


Re 


r 


norm 


S 

SA 

SST 

TLNS3D 

t 

tm 

URL 

U, V, w 


V 

K 


k y, z 


flap overhang, inches 
pressure, pounds /inch 2 
free-stream dynamic pressure. 


pM + 

2 


, pounds /inch 2 


Reynolds number based on the mean 
aerodynamic chord 

perpendicular distance above the model 
surface divided by the local chord of 
the wing in the cruise configuration 
semi-span wing area (i.e., reference 
area), inches 2 


Spalart-Allmaras turbulence model 

Menter's two-equation k-oa Shear-Stress- 

Transport turbulence model 

the Navier-Stokes code TLNS3D 

local wing thickness, inches 

turbulence model 

universal resource locator 

velocities in the physical coordinate x, y, 

and z directions, respectively, 

inches / second 


velocity magnitude, Vw 2 + v 2 + w 2 , 
inches / second 

free-stream velocity magnitude, 
V M x + U + K , inches / second 


physical (Cartesian) coordinates in the 
axial, spanwise, and vertical directions, 
respectively (origin at the intersection of 
the wing leading edge and the reflection 
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X 


X 


m/ 


Z 


m 


a 

6 

ti 

r\ norm 

T| + 

A 

P 

^ wall 

P 

Subscripts: 

c/4 

D 

f 

L 

le 

m 

norm 

s 


plane with the wing in the cruise 

configuration), inches 

axial distance from the leading edge of 

the main wing element divided by the 

local chord of the wing in the cruise 

configuration 

distances in the x and z directions. 


respectively from the (0, 0, 0) reference 
point to the pitching moment center, 
(see table 3), inches 
angle of attack, degrees 
deflection angle of leading-edge slat or 
trailing-edge flap, degrees 


wing span station, 


normal distance to the wall, inches 


law-of-the-wall coordinate, ^ norm ^^ > ' r — 

/d 


sweep angle, degrees 

slugs 

absolute viscosity, 

inch second 

local shear stress at wall, pounds /inch 2 
density, slugs /inch 3 


wing quarter chord 

drag 

flap 

lift 

wing leading edge 
pitching moment 
perpendicular to the model surface 
slat 
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V 


VISCOUS 


wall 

model surface 

wing 

wing 

00 

infinity, i.e., free stream conditions 


Multi-Element High-Lift Basics 

High-lift systems basically consist of trailing-edge devices that increase the 
turning angle of the flow over the wing, and thus increase lift. Slotted flaps and 
leading-edge slats additionally energize the boundary layer, and thus assist the 
flow in negotiating the greater turning angle without separating (see refs. 31 to 
34). Reference 31 briefly describes the milestones of their evolution from the 
aileron that Henry Farman was the first to use in 1908. 

Figure 1 depicts most of the well-known features of the flow about a well- 
configured multi-element wing at a low subsonic speed and a very high angle of 
attack. 2 Part (a) shows the overall flow; part (b) shows a detail of the flow in the 
vicinity of the leading-edge slat. Although very complex, much has been learned 
about its aerodynamic principles (see, for example, refs. 3, 7 to 18, 24, and 31 to 
36). Reference 35 gives a theoretical explanation of the effect of gaps, or slots, in 
multi-element airfoils for conditions at which the Reynolds numbers are 
sufficiently high to ensure that the boundary layers are turbulent by the time 
adverse pressure gradients appear. It lists the following primary effects for a 
properly designed and spaced combination of a leading-edge slat and other 
airfoil elements: 


2 Figure 1 presents a chord wise cut of a Navier-Stokes solution for a multi- 
element wing at a low subsonic speed and an angle of attack slightly greater than 
the angle for maximum lift. The lines illustrating the pattern of the flow are 
computed streamlines that have been numerically constrained to the plane 
presented. Throughout the paper, lines of this type will be referred to as 
"constrained streamlines". 
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(1) Slat effect: The circulation on the leading-edge slat opposes the 
circulation on the downstream elements. It hence reduces the negative 
pressure peaks on the downstream elements and in doing so relieves the 
pressure rises required of the boundary layer of the downstream 
elements. 

(2) Circulation effect: At the same time, the circulation on the downstream 
element increases the velocity at the trailing edge of the slat, which is 
inclined to the wing chord. This, therefore, increases the circulation on 
the slat. 

(3) Dumping effect: The trailing edge of a forward element, such as a slat or 
a main wing element, is in a region of higher velocity and hence lower 
pressure than free stream. This relieves the pressure rise required of the 
boundary layer of the forward element and thus reduces the chance of 
separation. 

(4) Off-the-surface pressure recovery: The boundary layer emanating from 
a forward element is dumped at higher velocities than free stream and is 
dumped away from the wall of the downstream element. Deceleration 
out of contact with a wall is much more efficient than deceleration in 
contact with a wall. 

(5) Fresh-boundary-layer effect: The boundary layer of each wing element 
starts out fresh at the leading edge. Hence the boundary layer on each 
element is relatively thin and can take a more adverse pressure gradient 
without separating than if it were thick. 3 

As stated previously, these effects can result in very complex flows, 
especially at high angles of attack (see refs. 3, 7, and 36; figure 1). First, the 

3 It should be reiterated that in reference 35, the boundary layers are considered 
to be turbulent by the time adverse pressure gradients appear. When a boundary 
layer starts fresh, it starts at a low Reynolds number and is thin. In discussing 
effect no. 5, reference 35 uses Stratford's turbulent separation criteria to show 
that a thin boundary layer can sustain a more adverse pressure gradient than a 
thick one. 
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constrained streamlines indicate that strong streamline curvature exists at the 
leading edge of each element of the wing. At the leading edge of the slat, the 
yellow-orange-red color in the figure indicates that transonic flow can develop. 
Each element of the wing develops its own boundary layer under the influence 
of the local stream; and, on the top surface of the wing the boundary layers grow, 
shed, overlap, diffuse, and finally merge. 

The figure also depicts separation of the flow from the surface in the slat 
cove, the wing cove, and at the trailing edge of the slat. (At lower angles of 
attack, the flow on the trailing edge of the flap is usually separated and the flow 
on the trailing edge of the slat is attached (ref. 27)). Off-body separation, or flow 
reversal, above and downstream of the airfoil is indicated, a phenomena that 
occurs at very high angles of attack. Note that a layer of attached flow remains 
underneath the off-body separation. These complexities plus boundary layer 
transition from laminar to turbulent flow, the possibility of relaminarization of 
the flow in regions of rapidly expanding flow, and attachment-line transition 
effects (see ref. 2), pose a strong challenge for a CFD code to accurately model. 


High-Lift Configuration 

The configuration selected for the joint computational and experimental 
project referred to in the introduction was dictated by a desire for one that was 
relatively simple, yet one that still produced the complex high-lift flows of 
current transport aircraft (reference 13). These requirements plus experimental 
considerations led to a semi-span, low-aspect-ratio, swept, trapezoidal wing 
mounted on a body pod. The wing was composed of three elements: the main- 
wing element, a leading-edge slat, and a trailing-edge flap. It was untwisted and 
had no dihedral. It has become known as the NASA Trapezoidal Wing. 

Figure 2 depicts the computational rendition of the model together with the 
grid topology and the computational domain for this study. The wind tunnel 
walls were not modeled. Figure 3, photographs of the experimental model 
installed in the Langley 14x22 Foot Wind Tunnel (ref. 13), shows the 
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configuration in greater detail. Part (a) of the figure shows a top view of the 
complete wing-body, whereas part (b) shows a bottom view at the wing-body 
juncture. In part (b), some of the brackets that attach the slat and flap to the main 
wing element can be seen. The brackets were not simulated in the numerical 
computations for this paper. 

Figure 4 presents sketches of a typical airfoil section of the wing with the 
slat and flap in both stowed (cruise) and deployed (high-lift) positions. It also 
lists dimensions of the wing 4 . The high-lift configuration (i.e., with the slat and 
flap deployed) was the only one considered in this paper; the cruise 
configuration is used only to report the dimensions of the wing. For the present 
investigation, the leading-edge slat was deployed at 30.0°, and the full-span 
trailing-edge flap was deployed at 25.0°. These settings are typical of a landing 
configuration. The width of the leading-edge slat is constant. The ratio of its 
chord (measured normal to its leading edge) to the mean aerodynamic chord of 
the complete wing is 0.126. The ratio of the flap's local chord to the local chord of 
the complete wing is constant at 0.3 throughout the span of the wing. The 
positions of the slat and flap relative to the main wing element for the present 
study are given in table 1. Sketches of airfoils at selected span stations for this 
configuration are given in figure 5. The sketches illustrate that since the chord of 
the leading-edge slat is constant, the slat comprises a larger percentage of the 
chord of the wing at the wing tip than it does at the wing root. 

Reference 22 gives more details of the configuration and the experimental 
model. A complete description of the model is given at the NASA Trapezoidal 
Wing Web Site (see the appendix for the URL of the site). 


4 Since the wing was designed so that the slat and flap could be mounted 
at various settings, the dimensions of the wing are given with the slat and flap in 
the cruise position, i.e., with the leading-edge slat and the trailing-edge flap 
stowed. 
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Numerical Procedure 


Computational Domain and Grid 

As stated previously, figure 2 illustrates the computational domain and the 
general topology of the grid. The directions of the physical coordinates of the 
grid are included. The domain and grid were configured as a collection of C- 
grids wrapped around the body pod and the elements of the wing (i.e., C-grids 
in the x-z plane). In the plane perpendicular to the x axis (i.e., the y-z plane) the 
grid assumed a pseudo-H topology. Symmetry was assumed about the y = 0 
plane. The origin of the coordinate system was defined as the intersection of the 
leading edge of the wing in the cruise configuration with the y = 0 plane. 

The domain extended 6.44 mean-aerodynamic-chord (MAC) lengths in 
front of the origin, 8.55 MAC lengths behind it, 6.35 MAC lengths above and 
below it, and 5.17 MAC lengths (2.5 wing semi-span lengths) normal to the 
reflection plane. For comparison, the body pod is 3.00 MAC lengths long. In 
retrospect, the extent of this domain may be too small for the 3-D high-lift flow 
field being investigated here. Grids used for more recent NASA Trapezoidal 
Wing computations extended about 100 MAC lengths from the wing. The total 
domain is subdivided into 46 grid blocks to facilitate a body-fitted grid around 
this configuration. Block-to-block connections consisted of both point-to-point 
matching and patched block surfaces. 

Three different grids were investigated: one with a total of 12.0xl0 6 nodes 
(fine grid), one with a total of 6.0xl0 6 nodes (medium grid), and another one 

with a total of 1.5xl0 6 nodes (coarse grid). Each grid has the same fundamental 

topology. The grids are body fitted (i.e., grid lines coincide with the model 
surface and other boundaries) to facilitate implementation of the boundary 
conditions. Grid lines are clustered in the vicinity of the model, normal to its 
surface, in the vicinity of each element of the wing, and near the wing tip. For the 
medium and fine grids, the number of grid planes in each computational 
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direction of each block supports three levels of multi-grid. The coarse (1.5xl0 6 - 
point) grid was created from the fine grid by dropping every other grid plane in 
each of the three computational directions. Therefore, it can support only two 
multi-grid levels. 

The fine grid was created to address problems (to be discussed 
subsequently) with the 6.0xl0 6 -point grid by judiciously refining it in regions of 
high flow gradients. The mesh density in the circumferential direction (i.e., 
around the airfoil) for each element of the wing was increased in order to better 
resolve the gradients in the nominal direction of the flow. The density in the 
spanwise direction near the wing tips was increased to better capture the wing 
tip vortices. The density normal to the model surfaces was also increased with 
grid points clustered much closer to the surfaces to better resolve the boundary 
layers. As a result, for this fine grid, the first grid point off the surface was 
approximately 0.00002 inches normal to the surface, which resulted in a 
calculated law-of-the-wall coordinate, r\ + , of approximately 0.5 over most of the 
surface of the wing with attached flow. The fine grid also had significantly lower 
stretching factors than the medium grid. It was eventually chosen as the baseline 
grid for the present study. It should be noted that the "fine" grid used for these 
older computations (circa 2004) would be considered "coarse" by the standards 
of today (circa 2014). 

Table 2 and figs. 6 and 7 present some details of the fine grid. The grid is 
composed of 46 grid blocks containing a total of 11,971,510 grid points and 
11,188,288 cells (See table 2.). Fig. 6 depicts the model surface grid with every 
other grid line deleted for clarity. Fig. 7 (a), (b), and (c) respectively illustrate the 
overall grid topology in a vertical cut through the wing, and depict the details of 
the mesh in the coves of the leading-edge slat and the main wing element. Parts 
(b) and (c) of the figure show that the small bases at the trailing edges of the 
leading-edge slat and the main wing element were computationally simulated. 
The base at the end of the trailing-edge flap was not simulated since it was much 
smaller than the other two bases. 
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Numerical Algorithm 

Numerical calculations for the present investigation were made with the 
three-dimensional thin-layer Navier-Stokes code TLNS3D (refs. 24 to 28). As a 
generalized thin-layer code, TLNS3D retains the diffusion terms along all 
principal directions, neglecting only the cross-diffusion terms. To advance the 
solution in time to a steady-state solution, TLNS3D utilizes a semi-discrete cell- 
centered finite-volume scheme based on a Runge-Kutta time stepping scheme. 
The code employs grid sequencing and multi-grid. Linear, fourth-difference and 
nonlinear, second-difference artificial dissipation terms are added to respectively 
suppress odd-even decoupling and oscillations in the vicinity of shock waves 
and stagnation points. Implicit residual smoothing is used to increase the 
allowable time step. TLNS3D models the effects of turbulence through an eddy- 
viscosity computed by either the one-equation Spalart-Allmaras or by Menter's 
two-equation Shear-Stress Transport (SST) turbulence model. 

Turbulence Models 

For all the conditions of this investigation, turbulent dissipation was 
modeled by the Spalart-Allmaras one-equation turbulence model (ref. 29). That 
model was chosen because it is well established, gives very good results over a 
wide range of flow conditions, and is available in most modern Navier-Stokes 
codes. For one angle of attack, a = 29.877°, Menter's shear-stress-transport (SST) 
turbulence model (ref. 30) was also employed. It is a two-equation, blended k-s 
and k-co model. It also has a limiter that allows it to more accurately account for 
the transport of the principal shear stresses in adverse pressure gradients than 
the Spalart-Allmaras model. The Menter's SST model has been shown to produce 
more accurate results than the Spalart-Allmaras model for some configurations 
with reverse flow regions (see refs. 9 and 28). Thus it was thought that it would 
possibly improve the solution in regions of incipient separation, in separated 
regions, or in highly vortical regions. 
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Boundary Conditions 

Because the free stream is subsonic, Riemann invariants for a one- 
dimensional flow were used to calculate the flow variables at the inflow, outflow, 
and far-field boundaries of the computational domain (see fig. 2). Reflection 
boundary conditions were imposed on the y = 0 surface of the computational 
domain 5 . Finally, a no-slip boundary condition was imposed on all the surfaces 
of the model. 

Since a cell-centered finite-volume algorithm is used in the code, these 
boundary conditions are employed on a layer of ghost cells (auxiliary cells) that 
surround the physical boundaries of the faces of the grid blocks. At grid block 
interfaces, a second layer of ghost cells is also added to facilitate the evaluation of 
the fourth-difference artificial-dissipation terms across these boundaries. 

Numerical Solution Strategy 

In the current work, a 5-stage Runge-Kutta scheme was chosen as the 
primary solver in the TLNS3D code. Grid sequencing and multigrid with a V- 
type cycle were also utilized. All solutions were started with the entire flow field 
set at free stream conditions. Implicit residual smoothing was used for enhancing 
the stability bounds of the numerical scheme. 

The solutions were obtained in single precision mode. However, the 
interpolation coefficients that transfer information between adjacent patched grid 
block boundaries were computed in double precision. As a check on the accuracy 
of the single-precision-solver mode, one double-precision solution was obtained 
at an angle of attack of 29.877°. The two modes of precision gave nearly identical 
results for lift, total drag, and skin-friction drag. The double-precision mode, 
however, resulted in a more negative pitching moment that agreed better with 
experiment than the single-precision mode. 


5 A reflection boundary condition was imposed since the calculations are 
compared to experimental data that was obtained with semi-span experimental 
test techniques that minimized the effect of the tunnel wall boundary layer on 
the model, and closely simulated full-span flow (refs. 13 and 37). 
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Typically, in obtaining a solution with three levels of multigrid, 150 time 
steps on the coarse grid were run. This was followed by 150 multigrid cycles on 
the intermediate grid using two grid levels. The remaining multigrid cycles 
required for the solution to converge to a steady state were then run on the fine 
grid using three grid levels. 

Origin Computer Resources Required 

The numerical solutions were obtained on the NASA Ames Research 
Center's SGI Origin computer 6 . For the fine (12xl0 6 -point) grid, running in the 
multigrid mode with three grid levels required the following resources. 


Type resource 

Amount of resource required 

Single precision 

Double precision 

Number of CPUs 

20 

38 

Number of nodes 

10 

19 

Memory/ grid point (kb) 

0.355040 

0.753503 

CPU time/ grid point/ multigrid cycle (sec) 

0.000188 

0.000298 

Wall-clock time / grid point / multigrid cycle (sec) 

0.00000949 

0.00000945 


Hence a single-precision solution required 4.25 gigabytes of memory, and took 
between 564 and 1112 CPU hours depending on the number of multigrid cycles 
required for the solution to converge to a steady state. The equivalent wall-clock 
time took between 28.4 and 60 hours. However, since the solutions were not 
obtained in one continuous computer run, but were broken into several 
computer runs of approximately 150 multigrid cycles each, the typical solution 
took one to two weeks to obtain due to scheduling policies. 


6 At the time the numerical calculations were made, 2001 to 2004, the Origin 
computer was one of the most powerful at NASA Ames Research Center. 
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Experimental Data 


The experimental data presented in this paper is from the NASA 
Trapezoidal Wing website. It is listed in table 3. The data consists of that portion 
from the initial (1998) NASA Langley 14x22 Wind Tunnel investigation of the 
NASA Trapezoidal Wing. That investigation was originally reported in reference 
13. Since then, references 18, 22, and 38 have added valuable information about 
the model, the experiment, and the repeatability of the data. For completeness, 
the appendix very briefly describes the aspects of that investigation that are most 
helpful in assessing the comparison of the present numerical results with that 
experiment. 


Results 

Solutions were obtained with the Navier-Stokes code TLNS3D for the 
NASA Trapezoidal Wing, a high-lift, three-element wing similar to a transport 
wing at takeoff and landing. For this study, the leading-edge slat and trailing- 
edge flap of the wing were deployed at 30.0° and 25.0° respectively. The brackets 
that attach the leading-edge slat and the trailing-edge flap to the main wing 
element were not simulated. The study was conducted at a nominal free-stream 
Mach number of 0.20, and at eight angles of attack ranging from 10.004° to 
34.858°. The Reynolds number based on the mean aerodynamic chord of the 

wing was 4.3x10^. All calculations were made for free air. The resulting solutions 
were compared with experimental (i.e., wind-tunnel) data. 

Vetting the Solutions 

Reference 39 defines "code verification", "code validation", and "code 
calibration" as: 

Code verification: solving the equations correctly 
Code validation: solving the correct equations 
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Code calibration: adjustment of constants and parameters needed to fit 
experimental data 

Reference 37 states that, based on these definitions, code validation is an 
engineering problem whereas code verification is a mathematical problem. 

In the past, TLNS3D has given excellent results on a large variety of 
problems (see, for instance, refs. 8 and 24 to 28). Although this does not verify 
that the code solves every detail of the equations correctly in all cases, 
verification in the present paper will be limited to a discussion of the numerical 
convergence of the solutions. Validation of the code and turbulence models that 
were tested will inherently be assessed for the present problem as the numerical 
solutions are discussed and compared with the experimental data. Since the 
purpose of the present investigation was to evaluate the "off-the-shelf" 
performance of the code and turbulence models, there has been no attempt to 
calibrate the code by adjusting constants or parameters to fit the solutions to 
experimental data 

Numerical convergence. Tables 4 and 5 and figure 8 illustrate the 
numerical, or iterative, convergence for the calculations presented in this paper. 
The numerical solutions were deemed converged to engineering accuracy when 
the residual dropped at least two orders of magnitude and the values of the lift 
and drag coefficients changed one tenth of one percent or less in the last 100 
multigrid cycles. Table 4 lists all the cases run. It includes the numerical precision 
and the number of multigrid cycles run. Table 5 lists the lift, drag, and pitching 
moment coefficient values, and the percentage change in these coefficients in the 
last 100 multigrid cycles for all the cases. As the table shows, the changes of the 
lift, drag, and pitching moment coefficient values are equal to or less than 0.1% in 
almost all the cases. The exception was at the combination of 10.852° angle of 
attack and the medium grid. The pitching moment also failed to reach these 
criteria at 14.000° angle of attack. 

Figure 8 illustrates a problem that only showed up for the combination of 
the medium grid and an angle of attack of 14.000°. The figure compares the 
numerical convergence for both the medium and fine grids at this condition. The 
residual for the solution with the medium grid dropped almost four orders of 
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magnitude. However, the solution never actually settled out. 7 Figure 8 also 
shows that although the residual for the fine grid only dropped approximately 
two orders of magnitude, the lift and drag coefficients settled out to a relatively 
constant value. These coefficients were monitored for all the solutions for each of 
the three grids and all angles of attack, and they settled out to a relatively 
constant value in all cases except the one just discussed. 

Grid convergence. Reference 40 states that for a mathematically formal 
grid convergence study of a multidimensional problem with a single grid size 
measure, all the grids must belong to the same family, i.e., the aspect ratio of 
each successive grid has to remain constant. In addition, the solutions on each 
grid must be sufficiently converged, preferably with residuals reduced several 
orders of magnitude. 

Due to the computer resources available and time constraints for the 
present investigation, strictly meeting the previous requirements was not 
possible. However, solutions obtained on the three grids, the fine grid, the 
medium grid, and the coarse grid, were compared in an attempt to access grid 
convergence. See the section, "Computational Domain and Grid", for a detailed 
description of the grids. 

Table 6 presents the force-and-moment results for these solutions at four 
angles of attack: 10.004°, 10.852°, 14.000°, and 29.877°. Initially, the medium grid 
was created, and solutions were obtained at three angles of attack, 10.852°, 
14.000° and 29.877°. The first angle of attack, 10.852°, was mistakenly run instead 
of 10.004°. A comparison of these three solutions with the experimental data 
indicated that the grid should be improved. Thus the a=10.852° solution was not 
repeated at the correct angle of attack. It is, however, presented in table 6 to help 
assess grid convergence. 

Table 6 shows that the lift coefficients are very close for the fine and 
medium grids. In addition, the lift coefficient for the coarse grid is relatively 
close to lift coefficient of the other two grids at an angle of attack of 10.000°. It 
does not agree as closely at an angle attack of 29.877°, a result that might be 

7 This numerical iteration was terminated at 2000 multigrid cycles to conserve 
computer resources. 
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expected considering the strong flow gradients encountered at this condition and 
the coarseness of the lowest grid density. 

Based on the results presented in table 6, the 12xl0 6 node, or fine, grid was 
selected as the baseline grid for this study. The subsequent comparison of the 
numerical calculations with experimental data, and the discussion of the 
predicted flow characteristics of the problem are all based on solutions obtained 
with the 12.0xl0 6 node, or fine, grid. 

Forces and Moments 

Figure 9 and table 7 present a comparison of the computed and 
experimental lift, drag, and pitching moment coefficients. The values of the 
experimental data listed in table 7 are linearly interpolated (for either a or C L as 
the independent variable) from the experimental data given in table 3. 

Lift coefficient. Figure 9(a) presents the computed and experimental lift 
coefficients versus angle of attack. At 10.004° angle of attack, the lift coefficients 
agree within 0.99%. At 19.937°, they agree within 2.11%. See table 7 for the exact 
values. These values, in this linear range, yield slopes of the lift curves of 0.065 
for the computations, a slope that is within 4% of the experimental value. At 
max i mum lift, traditionally one of the most difficult parameters to determine 
accurately (ref. 30), the calculated lift coefficient is approximately 5.40% low. See 
the following values taken from table 7 at constant a. 


a t'-L, numerical 

10.004 

14.000 

19.937 

27.908 

28.120 

29.877 

32.993 


'L, experimental/ ' Vk ~'L, experimental 

-0.99% 

-1.28% 

- 2 . 11 % 

-3.41% 

-4.80% 

-4.80% 

-5.40% 


The predicted or calculated angle of attack for maximum lift is about 
32.993°. This is close to the experimental value of 33.7°. The lack of additional 
solutions at angles of attack near 33° prevents a more exact assessment of the 
computed angle of attack for maximum lift, or for stall. At angles of attack 
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greater than the angle for maximum lift, the computed lift drops as it should, but 
not as abruptly as the experimental data. Similar results were reported in ref. 11. 

Attempts were made to get better agreement between the numerical 
calculations and the experimental data at 29.877° angle of attack. As stated 
previously in the section on the numerical procedure, the solutions were 
obtained in single precision on the NASA Ames Origin computer complex. A 
double precision solution was tried but showed no improvement. Substituting 
Menter's SST turbulence model for the Spalart-Allmaras model also failed to 
improve the agreement. As will be discussed subsequently, the prediction of 
max i mum lift possibly may be improved by increasing the spanwise density of 
the grid near the wing tip. Another possibility may be by modeling the test 
section of the wind tunnel (ref. 11). However, as previously stated, in the 
Langleyl4x22 Foot Tunnel, the wall interference is believed to be relatively small 
and the corrections to free air accurate. Even though the computed maximum lift 
is low, it is at least a reasonable estimate of the experimental value. 

The present high-lift problem has recently been the subject of many other 
investigations. Reference 41, for example, reports a very careful study of grid 
resolution for it, and of modeling the brackets that attach the slat and flap to the 
main wing element. By strategically refining the grid, it greatly improved the 
agreement between experiment and computations at high angles of attack. This 
was confirmed in reference 42 which included even finer grids. Including the 
attachment brackets decreased the computed C L and increased C D . Also see 
references 18, and 20 to 22 that summarize the results of the First AIAA CFD 
High-Lift Prediction Workshop. Examples of additional studies of this problem 
are: effects of grid resolution, references 43 to 45; effects of modeling the 
brackets, references 46 and 47; effects of modeling transition, references 47 to 50; 
effects of thin-layer versus full Navier-Stokes solvers, reference 46; and the 
effects of hysteresis, reference 51. A complete list of publications directly related 
to the First High-Lift Workshop can be accessed at the URL for the workshop. 

See the appendix. 

Drag coefficient. Figure 9(b) presents drag coefficient versus angle of 
attack. Again, the major discrepancy between the computations and experiment 
is at the higher angles. However, the discrepancy at angles of attack of 19.937° 
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and below is now on the order of -6.11% to -7.97%, and about 10% at the higher 
angles. At maximum lift, it is -8.58%. As with lift coefficient, the drag coefficient 
is under predicted at all angles of attack. See the values below taken from table 7. 
It is well known that drag is more difficult to predict accurately than lift (see ref. 
6 ). 


a 

(^-D, numerical experimental) / experimental) 

10.004 

-6.11% 

14.000 

-7.08% 

19.937 

-7.97% 

27.908 

-10.06% 

28.120 

-10.09% 

29.877 

-9.57% 

32.993 

-8.58% 

Part (c) of figure 9 

presents a plot of drag coefficient as a function of lift 


coefficient. It shows, as does the following comparison from table 7 at constant 
values of C L , an overall improvement over the comparison versus constant a. 
Versus constant C L , the computed drag coefficient is under-predicted by -4.06% 
to -5.45% at the lower values of lift, and up to +8.59% near the maximum value of 

c L . 


c L 

(C -C 

V^"D, numerical ^ D, experimental 

1.806 

-4.55% 

2.078 

-5.45% 

2.447 

-4.06% 

2.790 

0.62% 

2.796 

0.72% 

2.845 

2.32% 

2.867 

8.59% 


)/c 


D, experimental 


) 


That the numerical computations were made with the boundary layers fully 
turbulent leads one to question how this might effect the agreement between the 
computations and experiment. The fact that skin friction drag is such a minor 
part of the total drag (see table 5), suggests that skin friction, in this case, is of 
secondary importance. More importantly, modeling boundary layer transition 
would insure that the boundary layers on each element were in the correct state 
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as they mutually interacted. Transition also has a major influence on reducing 
flap separation (see refs 43, and 47 to 49) which affects the overall solution, 
including the drag. 

In addition to not precisely modeling the state of the boundary layer over 
the entire surface, the brackets that attach the slat and flap to the main wing 
element were not modeled for the computations. One would expect that 
including the brackets in the computations would increase the computed drag 
and thus the agreement with experiment. Again, reference 41 indicates that 
including the brackets in the computations significantly increased the drag but 
decreased the lift except at the very low angles of attack (see also refs. 20, 46, and 
47). 

Pitching moment coefficient. The calculated pitching moment coefficient, 
fig. 9(d), follows the correct trends with angle of attack. Similar to the lift 
coefficient, it agrees best with experiment for the angles of attack that fall within 
the linear range of the lift curve ( a = 19.937° or less). At the higher angles of 
attack, the pitching moment coefficient fails to compare well with the data (see 
table 7). This would be expected since both the lift and drag coefficients also fail 
to match the experimental data well at these angles. Reference 46 indicates that 
including the full viscous flux operator in Navier-Stokes solutions can 
significantly increase lift and reduce pitching moment. 


a v^-rn, numerical 

10.004 

14.000 

19.937 

27.908 

28.120 

29.877 

32.993 


'm, experimental/ / ^"m, experimental. 

-4.31% 

-5.24% 

-7.22% 

-13.20% 

-13.24% 

-13.16% 

-15.50% 


Force-and-moment summary . The present numerical study gave results 
similar to many of the entries submitted to the first AIAA CFD High-Lift 
Prediction Workshop, AIAA HiLiftPW-1 (ref. 20). The numerically computed lift, 
drag, and pitching moment of the present study exhibit the correct trends with 
angle of attack. The computed lift agrees relatively well with experiment in the 
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linear range (a = 19.937° or less) of the lift curve. In addition, the computations 
give reasonable estimates of the maximum lift and of the angle of attack at 
max i mum lift. These are all very important parameters in aircraft design that 
affect the efficiency, range, takeoff characteristics, and maximum payload of the 
airplane (refs. 2 and 30). Thus, the present study adds to the volume of work that 
helps confirm, and sheds insight into, the basic technology for calculating 
complex high-lift flows. 

Based on this basic technology, and with the help of more powerful 
computers, significant progress in predicting wind tunnel results has been made 
since the present calculations were completed. Again, see reference 20. 
Improvements due to modeling transition, improved grid resolution, modeling 
the brackets that attach the slat and flap to the main wing element, and including 
a full viscous flux operator (i.e., not thin-layer) in Navier-Stokes solvers have 
already been noted in the previous sections. These newer results imply that 
much of the disagreement between calculations and experiment for high-lift 
flows that has traditionally been attributed to turbulence modeling can be 
overcome by using sufficiently fine grids and by including these other details. 

Surface Pressures 

Figure 10 presents a comparison of the computed and experimental 
pressure coefficient distributions at several span stations of the wing. The 
computations presented in this figure were all made with the basic grid and 
turbulence model for the investigation (the fine grid and the Spalart-Allmaras 
turbulence model). The sketches in the figure depict the locations of the 
numerical and experimental pressures. The pressure distributions are presented 
at each angle of attack for which solutions were obtained (i.e., angles of attack 
ranging from 10.004° to 34.858°). 

Figure 10 illustrates that, over most of the span of the wing, there is 
generally very good agreement between the computed pressure distributions 
and experiment. For the leading-edge slat and the main wing element, the 
agreement is excellent at all angles of attack through the angle for maximum lift, 
32.993°. At the outermost span station of r\ = 0.95, the agreement between the 
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numerical computations and experiment is not as good. This exception will be 
discussed subsequently. 

Figure 10 further shows that the numerical computations correctly predict 
several trends that are indicated by the experimental data. First, as angle of 
attack is increased up to the angle for maximum lift, the slat and the main wing 
element become more heavily loaded, resulting in the negative pressure peaks on 
these elements increasing dramatically. In contrast, the negative pressure peak 
on the trailing-edge flap at first remains relatively constant with an increase in 
angle of attack, and then decreases at the higher angles. Also, as the angle of 
attack increases beyond the angle for maximum lift and the main element 
unloads, the unloading effect starts at the wing tip and progresses inboard. The 
calculations correctly capture this trend. However, the calculations predict that 
the trend starts at the slightly lower angle of attack of -34.0° vs. 35.5° for the 
experiment. 

Flow-Field Observations 

Conceptual illustration. Figure 11 conceptually illustrates simulated oil 
flows on the model surface as well as details of the flow in a vertical plane 
cutting the wing with the model at a very high angle of attack (~ 35°). Similar 
details of the computed flow immediately above the model surface and in this 
plane are respectively presented for each of the angles of attack investigated. See 
figures 12 to 19. These figures help illustrate the interaction of the boundary 
layers and wakes emanating from each of the three elements of the wing. 

Separation and correlation with surface pressures. Except at the outermost 
span station of q = 0.95, neither the computed nor experimental pressure 
coefficient distributions (figure 10) indicate much evidence of separation. 8 In 
contrast to the pressure coefficients, at a = 10.004°, the computed streamlines that 

8 The separation indicated by the experimental pressures at span stations of 0.50 
and 0.65 and at angles of attack above 28.120° is probably strongly influenced by 
the slat and flap attachment brackets that are located at q=0.47 and r|=0.61 
respectively. 
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are constrained very close to the surface of the trailing-edge flap show a small 
separated region at the trailing edge of the flap. The computations further predict 
that this separated region grows with angle of attack, reaches a maximum at 
about a = 19.937°, diminishes, and disappears altogether at a = 34.858°. Figure 20 
summarizes this effect. The computations also predict that as the flow on the 
flaps evolves from separated to attached flow, the flow on the leading-edge slat 
evolves from attached to separated flow. These trends are typical of multi- 
element high-lift systems that have been optimized for maximum lift (refs. 5 and 
27). Further examination of figs. 12 to 19 shows that this separation on the flap is 
limited to a very thin region close to the surface. (See the Mach number contours, 
the constrained streamline patterns, velocity vectors, and boundary layer profiles 
presented in parts (b), (c), and (d) of these figures.) Because the separation is so 
thin, it is not readily evident by examining the pressure distributions. At the very 
high angles of attack there is massive off-body separation that will be discussed 
subsequently. 

Shear layers and velocity profiles. The computed velocity vectors and 
velocity profiles close to the model surface also illustrate the various shear layers 
trailing from the different elements of the wing. (At the time when the present 
calculations were made, no experimental boundary layer profiles were available 
for comparison. Experimental velocity profiles for this configuration are now 
published in reference 38.) Below an angle of attack of 19.937°, the velocity 
vectors and profiles on the main wing element, parts (c) and (d) respectively of 
figures 12 to 19, do not show a significant velocity defect from the leading-edge 
slat. (The imaginary rakes depicted on the figure are to only show the location of 
the velocity profile. No physical rake was simulated for this numerical 
investigation.) However at a = 19.937°, the boundary layer on the main wing 
element (i.e., at x = 0.50) begins to show a defect that is probably the shear layer 
from the leading-edge slat. At angles of attack of 28.120° to 32.993°, the shear 
layer emanating from the slat is quite evident in the profiles on the main wing 
element. 

In contrast to the profiles on the main wing element, the boundary layer 
profiles on the trailing-edge flap clearly show the shear layer emanating from the 
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main wing element at all angles of attack. In addition, at angles of attack of 
28.120° to 32.993°, the shear layer from the slat, although it has largely dissipated, 
still appears to be evident in the profiles on the flap (at x bar = 0.85). Above a = 
32.993°, the profiles at both axial stations illustrate off-body separation with a 
thin underlying layer of attached flow. 

An examination of the details of the shear and boundary layers associated 
with each element of the wing suggests that each layer must be well resolved to 
accurately predict the aerodynamic force and moments. Reference 39 and 42 
confirm this speculation. Reference 39 showed that as angle of attack increases, 
properly resolving the boundary and shear layers emanating from the forward 
elements becomes increasingly important, especially at angles of attack near stall. 
Reference 48, for example, further showed that transition location and including 
rotation and curvature effects in turbulence modeling significantly affected the 
accuracy of predicting the forces and moments. References 43, 47, and 49 to 50 
also illustrate that transition has a major influence on flap separation. The 
complexity implies that automatic grid refinement could be very useful in 
adequately resolving the flow. 

Stagnation points. The computations show an interesting change in the 
stagnation line as angle of attack increases. As would be expected, and as figs. 12 
through 19 show, the stagnation line on the leading-edge slat moves aft as angle 
of attack is increased. It moves from approximately the leading edge of the slat to 
the corner of the slat cove where it appears to set up a saddle point just off the 
surface. Also, as would be expected, the stagnation line of the main wing element 
moves aft with increasing angle of attack. However, for the main wing element, 
the location of the stagnation points near the wing root and tip remain relatively 
constant with increasing angle of attack. This results in a rearward bowing of the 
stagnation line in the middle of the wing. In contrast, the stagnation line on the 
leading-edge slat remains relatively straight. Figure 21 shows computed oil flows 
(actually streamlines constrained immediately above the model surface) for the 
bottom view of the juncture of the leading edge of the wing with the fuselage. 
Unlike the slat and main wing element, the stagnation line for the trailing-edge 
flap moves very little with changes in angles of attack (refer to figs. 12 to 19). 
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Off-body separation. Figures 17 to 19 show that for the computations, at 
the angle of attack of maximum lift, a = 32.993°, the flow begins to separate at the 
outboard section of the leading-edge slat. In addition, the flow develops a 
massive low velocity region off the body that is centered just downstream and 
above the trailing-edge flap. As the angle of attack increases, the separation on 
the leading-edge slat progresses inboard. Simultaneously, the low velocity region 
moves forward, off-body separation develops, and the accompanying vortices 
develop in the huge separated region. Underneath this huge separated region, a 
thin layer of attached flow remains next to the upper surfaces of the main wing 
element and the trailing-edge flap. See the velocity vectors and the boundary 
layer profiles, parts (c) and (d) of the figures. 

Wing-tip vortices. Figure 22, which presents plots of the wing vortices at 
several angles of attack, also shows these trends of separation with angle of 
attack. At a = 10.004°, the wing tip vortex and the vortex sheet behind the 
trailing edge of the wing are well defined. At a = 29.877°, the tip vortex and 
vortex sheet begin to dissipate. At a = 32.993°, vortex dissipation increases, and 
separation at the tip of the leading-edge slat can be detected. (Note the 
streamline that originates on the slat and the flow pattern depicted in fig. 17(a).) 
At a = 34.858°, massive separation occurs and the tip vortex is well dissipated. 

As one would expect, this complex flow at the wing tip is extremely difficult to 
predict. Other research as reported in at the First AIAA HiLift Prediction 
Workshop (ref. 40) confirms this difficulty. As a matter of interest, a vortex 
forming at the juncture of the wing and body, another well-known phenomena 
(ref. 27), can be seen to develop and become well defined as angle of attack 
increases. 

These phenomena are reflected in the computed pressure coefficients 
(compare fig. 10 parts (e), (f), (g), and (h)). As angle of attack increases beyond 
the angle for maximum lift ( a ~ 33.0°), the peak expansions progressively become 
less severe. Correspondingly, the recompressions that occur over the latter part 
of the airfoil progressively become flatter and less negative. These effects start at 
the wing tip and progress inboard as angle of attack increases. The experimental 
data tend to exhibit the same pattern, although its onset is delayed. In addition. 
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near the tip of the wing, the levels of the pressure coefficients on the leading- 
edge slat also drop and flatten out, thus indicating separation on the slat. Again, 
the fact that the numerical procedure correctly predicts these data and known 
trends (refs. 5 and 24) reflects very positively on its ability to handle very 
complex flows of this nature. 

At the outermost span station, r\ = 0.95, the computed pressure coefficients 
on the trailing-edge flap are less negative than experiment (see fig. 10). This 
discrepancy between computations and experiment increases with increasing 
angle of attack. The increase in the discrepancy eventually works its way 
forward onto the main wing element and even onto the leading-edge slat at a = 
33.980° and 34.858°. Since this span station is in the vicinity of the wing tip vortex 
(see fig. 23), the discrepancy is probably due to the failure of the computations to 
accurately model the wing- tip vortex. 

As this low velocity region is moving forward with increasing angle of 
attack, the influence of the wing tip vortex and off-body separation is 
progressing from the wing tip inboard. Therefore, it is possible that these 
phenomena, which are simultaneously changing with angle of attack, were 
changing at slightly different rates or conditions for the computed and 
experimental flows. Under these circumstances, precisely matching the 
experimental results with numerical calculations would almost be fortuitous. 

Other possible sources of error include inadequate grid density (see refs. 

20, 41, and 42). Since the locations of the tip vortex and the shear layers trailing 
from the individual wing elements change with angle of attack (see figs. 12 to 19, 
and 22), the areas of the grid that are refined to capture these features should 
track them as angle of attack changes. Modeling the leading-edge slat and 
trailing-edge flap brackets, and modeling boundary layer transition may also 
help in pinpointing the source of these differences near the wing tip. Later 
research (ref. 46) suggests that Navier-Stokes algorithms that include the cross 
derivative diffusion terms, terms not included in the present thin-layer 
computations, can significantly reduce separation at the wing tip and improve 
agreement with experiment. As always, one wonders if the turbulence model is 
performing adequately in this particularly complex flow. 
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Effect of the Turbulence Model 

A very limited investigation was conducted on the effect of the turbulence 
model on the solution. At the angle of attack of 29.877°, in addition to the 
Spalart-Allmaras one-equation turbulence model, Menter's two-equation SST 
turbulence model was investigated. Menter's SST model is a blended k-e k-to 
model. It utilizes a limiter that allows it to more accurately account for the 
transport of the principal shear stresses in adverse pressure gradients than the 
Spalart-Allmaras model. Hence, the author thought that Menter's model might 
capture the physics better in regions of pending separation, separation, or maybe 
even strong vortical flow than the Spalart-Allmaras model. However, as table 4 
and figure 9 show, Menter's model yielded a two percent lower lift coefficient 
value than the Spalart-Allmaras yielded, and therefore did not improve the 
correlation of the computations with the experimental data. 

Figure 23 presents the pressure distributions for the two turbulence models 
at span stations of r\ = 0.70 and r\ = 0.95. Little difference exists between the 
distributions for the two turbulence models except for a slight difference on the 
trailing-edge flap. However, these seemingly insignificant differences in 
pressures can integrate into significant differences in forces as they have in the 
present case. See figure 9 and tables 4 and 5 (cases 7 and 7-SST). 

Figure 24 presents the effect of the turbulence model on details of the 
computed flow. There is relatively little difference in the constrained streamlines 
immediately above the surface (part (a)) between the two turbulence models 
except in a small region near the trailing edge of the trailing-edge flap at the 
wing tip. This might lead one to conclude that the turbulence model makes 
relatively little difference in the solution. However, an examination of the large 
region of low velocity flow downstream of the trailing-edge flap shows that the 
constrained strea mlin es and velocity vectors for the two models (parts (b) & (c)) 
exhibit a significant difference in the solutions in this region. 

This difference is further illustrated by the velocity profiles presented in 
figure 25. At the span station r| = 0.70, which is typical of all inboard span 
stations, there is little difference between the velocity profiles predicted by the 
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two turbulence models on the main wing element. The profile immediately next 
to the surface of the trailing-edge flap is approximately the same for the two 
turbulence models also. However, the shear layer emanating from the main wing 
element is slightly thicker and at a lower velocity for Menter's SST model than 
for the Spalart-Allmaras model. 

At the outermost span station, r\ = 0.95, the profiles indicate a huge low- 
velocity region over both the main wing element and the trailing-edge flap. As 
mentioned before, there is an underlying high-velocity layer of attached flow 
next to the surface. In this case, Menter's SST turbulence model gave a higher 
velocity than the Spalart-Allmaras model. In addition, the underlying high- 
velocity layer of attached flow on the flap is thicker and has a decidedly different 
profile. Perhaps these differences are an indication of a significantly different 
wing-tip vortex generated by the two models. See figure 26, which illustrates the 
effect of the turbulence models on the wing vortices. 

Figure 26 shows that the Menter two-equation SST turbulence model gives 
a slightly better developed wing tip vortex than the Spalart-Allmaras one- 
equation turbulence model. Also the vortex sheet immediately behind the wing 
appears to be stronger for Menter's model. These results indicate that the 
turbulence model can have a significant effect on the solution of three- 
dimensional multi-element high-lift flows. Reference 48, an investigation of 
transition and advanced turbulence modeling for this problem, tends to confirm 
this hypothesis. It shows that including transition and rotation/ curvature effects 
usually improved correlation with experiment. Again, see reference 46 which 
indicates that including the full viscous flux operator in Navier-Stokes solutions 
can significantly increase lift and reduce pitching moment. 


Concluding Remarks 


Solutions were obtained with the Navier-Stokes CFD code TLNS3D to 
predict the flow about the NASA Trapezoidal Wing, a high-lift wing composed 
of three elements: the main-wing element, a deployed leading-edge slat, and a 
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deployed trailing-edge flap. The brackets joining the wing elements were not 
modeled. Turbulence was modeled by the Sparlart-Allmaras one-equation 
turbulence model. For comparison, one solution at 29.877° angle of attack was 
repeated with Menter's two-equation k-co Shear-Stress-Transport (SST) 
turbulence model. All calculations were made with the assumption that the 
boundary layers were fully turbulent. The investigation was conducted at a free 
stream Mach number of 0.2, and at angles of attack ranging from 10.004° to 
34.858°. 

The numerical results predicted the correct trends of lift, drag, and pitching 
moment with angle of attack. In the linear range of the lift curve, a = 10.004° to 
19.937°, the computed and experimental lift coefficients agree well. At 10.004° 
angle of attack, the computations are 0.99% low, and at 19.937°, 2.11% low. This 
results in a predicted lift curve slope that is within 4% of the experimental value. 
At maximum lift, the calculated lift coefficient is approximately 5.40% low. The 
predicted or calculated angle of attack for maximum lift is about 33.0° versus an 
experimental angle of 33.7°. Versus angle of attack, the computed drag coefficient 
is between 6.11% and 7.97% lower than experiment in the linear range of the lift 
curve, and 8.58% lower at maximum lift. Versus C L , the drag coefficient is 4.55% 
to 4.06% lower in the linear range, and 8.59% higher at maximum lift. The 
predicted pitching moment, largely driven by the lift, is 4.31% to 7.22% lower 
experimentally at angles of attack within the linear range of the lift curve, and 
15.50% at the angle for maximum lift. As one might expect, at angles of attack 
beyond the angle for stall, all three of the predicted aerodynamic parameters fail 
to match the experimental data well. These force and moment results are similar 
to many other calculations for this problem. 

Over most of the span of the wing, there is generally very good agreement 
between the values of the computed and experimental pressure coefficients. For 
the leading-edge slat and the main wing element, the agreement is excellent at all 
angles of attack through the angle for maximum lift, 33.0°. At the outermost span 
station of q = 0.95, which is heavily influenced by the wing-tip vortex, the 
agreement between the computations and experiment is not as good. 
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The numerical computations correctly predict several trends and features of 
the flow that are typical of multi-element systems optimized for maximum lift. 
First, as angle of attack is increased up to the angle for maximum lift, the slat and 
the main wing element become more heavily loaded, with the negative pressure 
peaks on these elements increasing dramatically. In contrast, the negative 
pressure peak on the trailing-edge flap at first remains relatively constant with 
an increase in angle of attack, and then decreases at the higher angles. The 
computations predict that there is a separated region on the trailing-edge flap at 
low angles of attack. The separated region disappears at the higher angles. In 
contrast, the predicted flow on the slat is at first attached, and then separates as 
angle off attack is increased through the same range. Experimental data confirm 
these trends. The numerical procedure also predicts that at the angle of attack of 
maximum lift, the flow develops a massive low velocity region off the body. As 
the angle of attack increases, the low velocity region moves forward, off-body 
separation develops, and the accompanying vortices develop in the huge 
separated region. Underneath this huge separated region, a layer of attached 
flow remains next to the upper surfaces of the main wing element and the 
trailing-edge flap. These trends and features of the flow start at the wing tip and 
progress inboard with increasing angle of attack. 

At 29.877° angle of attack, Menter's SST turbulence model yielded a 2% 
lower lift than the Spalart-Allmaras model. It also gave a slightly better 
developed wing tip vortex than the Spalart-Allmaras model. In addition, the 
constrained streamlines and velocity vectors for the two models are significantly 
different in the large region of low velocity downstream of the trailing-edge flap. 

The overall good prediction of the lift parameters and surface pressures 
coupled with the correct prediction of many flow features and trends favorably 
reflects on the ability of present Navier-Stokes procedures to solve this complex 
high-lift problem. Since the present computations were made, many other strides 
(largely made possible by more powerful computers) have been made in this 
area. Using finer grids, a more detailed modeling of the configuration (such as 
modeling attachment brackets), including boundary layer transition, and 
including the cross-derivative diffusion terms (full Navier-Stokes) have all 
significantly improved the agreement between calculations and experiment. The 
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success of these later techniques suggests that their application can overcome 
much of the traditional discrepancy between computations and experiment, a 
discrepancy that has often been attributed to the turbulence model. However, 
progress is still necessary before numerical results can routinely replace wind 
tunnel data for high-lift flows. 
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Appendix 


Experimental Apparatus and Procedure 

The experimental data presented in the present paper is from the NASA 
Trapezoidal Wing Web Site. The force, moment, and pressure data presented is 
from the initial (1998) NASA Langley 14x22 Wind Tunnel portion of the 
investigation of the NASA Trapezoidal Wing. That investigation was originally 
reported in reference 13. Since then, references 18, 22, and 38 have added 
valuable information about the model, the experiment, and the repeatability of 
the experimental data. The "original" trap wing web site (i.e., the data for this 
paper) can be accessed from a link to the 1st Workshop (HiLiftPW-1) Website 
Home Page. The link is listed in the new HiLiftPW Home Page: 
http:/ / hiliftpw.larc.nasa.gov 

For completeness, this appendix briefly describes the aspects of that investigation 
that are most pertinent in helping assess the comparison of that experimental 
data with the present computational results. 

Wind Tunnel 

The Langley 14x22 Subsonic Tunnel is a closed-circuit, single-return, 
atmospheric wind tunnel with a maximum speed of about 338 ft/ sec. The 
max i mum unit Reynolds number is 2.1xl0 6 per foot. Low levels of test section 
turbulence are achieved through the combined effects of a fan diffuser grid, a 
flow-straightening honeycomb, four anti-turbulence screens, and a 9-to-l 
contraction ratio. To minimize test-section- wall interference on the flow over 
research models, the test section can be operated in a variety of configurations - 
closed, slotted, partially open, and open. The tunnel may also be configured with 
a floor boundary-layer removal system at the entrance to the test section. For the 
present investigation, the test section and slots were both closed. Reference 23 
gives a through description of the wind tunnel. 
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Model 

For the experimental portion of the combined project referred to in the 
introduction, a large model with a high Reynolds number and a thick boundary 
layer that could be measured with pneumatic instrumentation was highly 
desired (ref. 13). These requirements led to a semi-span, low-aspect-ratio, swept, 
trapezoidal wing mounted on a body pod. The wing was composed of three 
elements: the main-wing element, a leading-edge slat, and a trailing-edge flap. It 
was untwisted and had no dihedral. The wing of the experimental model had a 
semi span of 85.1 inches, and a mean aerodynamic chord of the wing of 39.6 
inches. Since the wing was designed so that the slat and flap could be mounted at 
various settings, the dimensions of the wing are given with the slat and flap in 
the cruise position, i.e., with the leading-edge slat and the trailing-edge flap 
stowed. Note that the high-lift configuration (i.e., with the slat and flap 
deployed) was the only one investigated in the present numerical study; the 
cruise configuration is used only to report the dimensions of the wing. 

The model was mounted from the floor of the wind tunnel as shown in 
figure 3. Mounted in this way, its wing tip extended to within approximately 6 
inches of the tunnel centerline. Figure 3 shows that the wing was outfitted with 
the full-span flap for this investigation (ref. 13). References 18, 22, and 38 list 
more information about the model, and a complete description of it is given at 
the NASA Trapezoidal Wing Web Site. The site contains files that provide the 
model surfaces in Initial Graphics Exchange Specification (IGES) format, and 
gives detailed instructions for transforming the slat and flap between their 
stowed and deployed positions. 

Instrumentation 

Force and moment data were obtained with a strain-gage balance that was 
accurate to within +/- 3.0 pounds normal force, + 1 - 2.2 pounds axial force, and 
+ / - 47.0 inch- pounds pitching moment (these are a 2-sigma level). The nominal 
free-stream dynamic pressure for this test was 58.3 pounds per square foot. 
Hence, based on the model reference area, etc. (figure 4), the typical balance 
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accuracies stated in terms of aerodynamic coefficients are approximately: + /- 
0.0024 in its normal direction, + /- 0.0017 in its axial direction, and + /- 0.00095 in 
pitching moment. The pressures were measured with electronic pressure 
scanning modules of 1.0, 2.5, 5.0, and 15.0 pounds per square inch full scale. Each 
pressure module was accurate to within 0.05 percent of full scale. 

Procedure 

The wind tunnel investigation (reference 13) was conducted at a free stream 
Mach number of 0.2, and at angles of attack ranging from -4.0° the 36.0°. To 
obtain a relatively thin tunnel-wall boundary layer in the vicinity of the model, 
the tunnel-floor boundary layer was removed at the entrance to the test section, 
(see ref. 37 for semi-span testing techniques.) In addition, the model was 
mounted with the y = 0 surface of the body pod slightly above the tunnel floor. 
Reference 37 states that for semi-span testing, a standoff distance of two times 
the wall boundary layer displacement thickness resulted in the closest simulation 
of full-span results. These techniques were followed to closely simulate a 
reflection condition at y = 0. 

The experimental force and moment data corrected for wind-tunnel wall 
interference by classical wall corrections are used in the present paper. Table 3 
lists the data. These Langley wind-tunnel data are considered largely 
interference free, and the corrections to free air are considered accurate (see refs. 
13, 23, and 38). The Trapezoidal Wing Web Site described previously gives 
detailed information for the experimental data and investigation. References 18, 
22, and 38 give a thorough description or the model, the investigation, and 
analysis of the data. They describe the repeatability of the data including the 
effect of hysteresis with increasing and decreasing angle of attack and the effect 
of tunnel entries. 
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Table 1. Settings for the Leading-edge slat and Trailing-edge flap 


Leading-edge slat 

6 S 

§ / ^local 

/ ^local 

30° 

0.015 

0.015 

Trailing-edge flap 

6 f 

§f/ ^local 

Of / C-local 

25° 

0.015 

0.005 


Notes: 


Slat: 

• Slat settings are specified in a vertical plane normal to the wing leading 
edge 

• Slat height (h s ) is measured relative to the surface hilite of the main wing 
element 

• Slat gap (g s ) is defined as the distance between the slat trailing edge and 
the surface of the main wing element 

• For a point on the trailing edge of the slat, the chord used to normalize the 
slat gap and height is the local stream-wise chord of the cruise planform 
passing through the identical point on the trailing edge of the slat with the 
slat and flap in the cruise (stowed) position 


• Flap setting are specified in a plane aligned with the free-stream (the x-z 
plane) 

• Flap gap (g f ) is specified as the distance from the trailing edge of the wing 
cove to the nearest point on the flap 

• Flap overlap (o f ) is specified as the steam-wise distance from the leading 
edge of the flap to the trailing edge of the cove 
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Table 2. Grid blocks, points, and cells for the fine grid. 


Grid 

Block 

Dimension 

No. 

Points 

No. 

Cells 

1 

J 

K 







1 

177 

9 

169 

269217 

236544 

2 

89 

65 

25 

144625 

135168 

3 

277 

85 

13 

306085 

278208 

4 

113 

17 

9 

17289 

14336 

5 

169 

21 

17 

60333 

53760 

6 

521 

13 

21 

142233 

124800 

7 

189 

29 

73 

400113 

379008 

8 

189 

29 

61 

334341 

315840 

9 

189 

29 

61 

334341 

315840 

10 

21 

45 

201 

189945 

176000 

1 1 

65 

53 

81 

279045 

266240 

12 

65 

53 

61 

210145 

199680 

13 

65 

53 

61 

210145 

199680 

14 

129 

29 

97 

362877 

344064 

15 

129 

29 

97 

362877 

344064 

16 

129 

29 

97 

362877 

344064 

17 

129 

29 

97 

362877 

344064 

18 

89 

21 

121 

226149 

211200 

19 

85 

21 

121 

215985 

201600 

20 

113 

17 

121 

232441 

215040 

21 

29 

69 

121 

242121 

228480 

22 

9 

25 

121 

27225 

23040 

23 

145 

29 

97 

407885 

387072 

24 

141 

29 

97 

396633 

376320 

25 

141 

29 

97 

396633 

376320 

26 

145 

29 

97 

407885 

387072 


Grid 

Dimension 



No. 

No. 

Block 

1 

J 

l< 

Points 

Cells 







27 

21 

57 

193 

231021 

215040 

28 

29 

21 

193 

117537 

107520 

29 

65 

29 

193 

363805 

344064 

30 

109 

29 

73 

230753 

217728 

31 

109 

29 

61 

192821 

181440 

32 

109 

29 

61 

192821 

181440 

33 

85 

21 

77 

137445 

127680 

34 

33 

33 

77 

83853 

77824 

35 

9 

49 

77 

33957 

29184 

36 

21 

81 

77 

130977 

121600 

37 

245 

21 

77 

396165 

370880 

38 

129 

29 

97 

362877 

344064 

39 

125 

29 

97 

351625 

333312 

40 

125 

29 

97 

351625 

333312 

41 

129 

29 

97 

362877 

344064 

42 

201 

25 

37 

185925 

172800 

43 

473 

13 

57 

350493 

317184 

44 

473 

13 

61 

375089 

339840 

45 

473 

13 

61 

375089 

339840 

46 

369 

73 

9 

242433 

211968 


TOTAL 


11971510 


1 1188288 
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Table 3. Wind-Tunnel Force-and-Moment Data 


a, 

degrees 

Moo 

_l 

U 

Cd 

Cm 

Xm / 

inches 

Z ID/ 

inches 








-3.936 

0.2038 

0.2838 

0.1139 

-0.2200 

34.342 

0.000 

-2.085 

0.2033 

0.4873 

0.1089 

-0.2665 

34.342 

0.000 

-0.096 

0.2021 

0.8223 

0.1182 

-0.3695 

34.342 

0.000 

2.084 

0.2022 

1.1944 

0.1442 

-0.4687 

34.342 

0.000 

4.009 

0.2027 

1.3578 

0.1741 

-0.4828 

34.342 

0.000 

6.029 

0.2018 

1.5261 

0.2075 

-0.4964 

34.342 

0.000 

8.026 

0.2013 

1.6804 

0.2420 

-0.5013 

34.342 

0.000 

10.004 

0.2019 

1.8240 

0.2783 

-0.5035 

34.342 

0.000 

12.020 

0.2016 

1.9701 

0.3187 

-0.5044 

34.342 

0.000 

13.969 

0.2017 

2.1027 

0.3614 

-0.5042 

34.342 

0.000 

15.928 

0.2022 

2.2416 

0.4071 

-0.5061 

34.342 

0.000 

17.970 

0.2015 

2.3744 

0.4549 

-0.5056 

34.342 

0.000 

18.916 

0.2025 

2.4347 

0.4777 

-0.5053 

34.342 

0.000 

19.937 

0.2029 

2.4997 

0.5030 

-0.5056 

34.342 

0.000 

20.925 

0.2016 

2.5603 

0.5275 

-0.5055 

34.342 

0.000 

21.999 

0.2023 

2.6243 

0.5543 

-0.5028 

34.342 

0.000 

22.914 

0.2029 

2.6789 

0.5769 

-0.4993 

34.342 

0.000 

24.014 

0.2031 

2.7460 

0.6062 

-0.4960 

34.342 

0.000 

25.005 

0.2025 

2.7986 

0.6313 

-0.4916 

34.342 

0.000 

25.977 

0.2028 

2.8460 

0.6546 

-0.4866 

34.342 

0.000 

26.941 

0.2025 

2.8885 

0.6788 

-0.4811 

34.342 

0.000 

27.908 

0.2022 

2.9310 

0.7017 

-0.4759 

34.342 

0.000 

28.970 

0.2023 

2.9652 

0.7225 

-0.4645 

34.342 

0.000 

29.877 

0.2027 

2.9884 

0.7401 

-0.4528 

34.342 

0.000 

30.910 

0.2029 

3.0096 

0.7585 

-0.4388 

34.342 

0.000 

31.319 

0.2028 

3.0161 

0.7654 

-0.4325 

34.342 

0.000 

31.922 

0.2034 

3.0242 

0.7767 

-0.4225 

34.342 

0.000 

32.587 

0.2028 

3.0284 

0.7855 

-0.4112 

34.342 

0.000 

32.993 

0.2029 

3.0306 

0.7917 

-0.4031 

34.342 

0.000 

33.665 

0.2031 

3.0283 

0.7982 

-0.3856 

34.342 

0.000 

33.894 

0.2029 

3.0259 

0.7994 

-0.3793 

34.342 

0.000 

34.452 

0.2033 

3.0143 

0.8054 

-0.3606 

34.342 

0.000 

34.858 

0.2034 

2.9884 

0.8112 

-0.3411 

34.342 

0.000 

35.401 

0.2060 

2.2707 

0.8210 

-0.2463 

34.342 

0.000 

35.996 

0.2054 

2.2482 

0.8302 

-0.2324 

34.342 

0.000 


Note: These data from the TrapWing web site have been corrected to free-air 
conditions. The values of C D on the site were updated on July 1, 2009. Previously, 
they were incorrect due to an error in applying wind-tunnel-boundary corrections. 
See reference 18. 
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Table 4. Computational Cases Run 


Case 

a, 

degrees 

Grid (No. 
of 

Points) 

Turbulence 

Model 

Numerical 

Precision 

(BITS) 

No. of 

Multigrid 

Cycles 

Run 

1 

10.004 

12.0M 

S-A 1 -equa. 

32 

900 

2 

10.050 

12.0M 

S-A 1 -equa. 

32 

900 

3 

14.000 

12.0M 

S-A 1-equa. 

32 

1135 

4 

19.937 

12.0M 

S-A 1 -equa. 

32 

1150 

5 

27.908 

12.0M 

S-A 1 -equa. 

32 

1425 

6 

28.120 

12.0M 

S-A 1 -equa. 

32 

1325 

7 

29.877 

12.0M 

S-A 1 -equa. 

32 

1325 

8 

32.993 

12.0M 

S-A 1 -equa. 

32 

1485 

9 

33.980 

12.0M 

S-A 1-equa. 

32 

1775 

10 

34.858 

12.0M 

S-A 1 -equa. 

32 

1625 

7-DP 

29.877 

12.0M 

S-A 1 -equa. 

64 

1305 

7-SST 

29.877 

12.0M 

SST 2-equa. 

32 

1600 

1-CG 

10.004 

1.5M 

S-A 1 -equa. 

32 

2000 

7-CG 

29.877 

1.5M 

S-A 1 -equa. 

32 

2000 

la-MG 

10.852 

6.0M 

S-A 1 -equa. 

64 

900 

3-MG 

14.000 

6.0M 

S-A 1 -equa. 

64 

2000 

7-MG 

29.877 

6.0M 

S-A 1-equa. 

64 

1700 


Note: 

S-A 1-equa. = Spalart-Allmaras 1 -equation turbulence model 
SST 2-equa. = Mentor's SST 2-equation turbulence model 
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Table 5. Calculated Forces and Moments with Numerical Convergence 


Case 

a, 

degrees 

Calculated Values of C L , C D v , C D , and C m & 

The Change in C L , C D , and C m in the Last 1 00 Multigrid Cycles 

C L 

AC l 

C D .v 

Cd 

AC d 

c m 

AC m 

1 

10.004 

1.806 

0.055% 

0.0124 

0.2613 

0.078% 

-0.4818 

0.124% 

2 

10.050 

1.810 

0.000% 

0.0124 

0.2622 

0.076% 

-0.4819 

0.083% 

3 

14.000 

2.078 

0.000% 

0.0126 

0.3365 

0.030% 

-0.4778 

0.042% 

4 

19.937 

2.447 

0.041% 

0.0127 

0.4629 

0.022% 

-0.4691 

0.064% 

5 

27.908 

2.790 

0.036% 

0.0124 

0.6311 

0.000% 

-0.4131 

0.097% 

6 

28.120 

2.796 

0.071% 

0.0124 

0.6346 

0.047% 

-0.4109 

0.024% 

7 

29.877 

2.845 

0.000% 

0.0122 

0.6693 

0.000% 

-0.3932 

0.127% 

8 

32.993 

2.867 

0.000% 

0.0116 

0.7238 

0.028% 

-0.3406 

0.088% 

9 

33.980 

2.664 

0.075% 

0.01 10 

0.7593 

0.026% 

-0.2669 

0.262% 

10 

34.858 

2.537 

0.1 18% 

0.0106 

0.7870 

0.063% 

-0.2359 

0.297% 

7-DP 

29.877 

2.845 

0.035% 

0.0122 

0.6693 

0.090% 

-0.3937 

0.025% 

7-SST 

29.877 

2.788 

0.036% 

0.01 18 

0.6483 

0.01 5% 

-0.3733 

0.000% 

1-CG 

10.004 

1.730 

0.058% 

0.0131 

0.2666 

0.000% 

-0.4668 

-0.021% 

7-CG 

29.877 

2.542 

0.000% 

0.01 18 

0.6457 

0.002% 

-0.3362 

0.000% 

la-MG 

10.852 

1.803 

0.222% 

0.0127 

0.2687 

0.186% 

data lost 

data lost 

3-MG 

14.000 

2.038 

0.098% 

0.0130 

0.3399 

0.088% 

data lost 

data lost 

7-MG 

29.877 

2.845 

0.035% 

0.0122 

0.681 1 

0.029% 

data lost 

data lost 
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Table 6. Grid Convergence Study 


Grid, no. 
of cells 

Angle of attack, a 

10.004° 

10.852° 

14.000° 

29.877° 

C L 

Cd 

C L 

Cd 

C L 

Cd 

C L 

Cd 

1.5x10 6 

1.730 

0.2666 

no data 

no data 

no data 

no data 

2.542 

0.6457 

6.0x1 0 6 

no data 

no data 

1.803 

0.2687 

2.038 

0.3399 

2.845 

0.6811 

1 2.0x1 0 6 

1.806 

0.2613 

no data 

no data 

2.078 

0.3365 

2.845 

0.6693 


Note: Due to the resources available, a mathematically formal grid study could not be 
conducted. These cases were run in order to make an engineering assessment of grid 
convergence. Refer to the section “Grid convergence” under “Results” in the text. 
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Table 7. Comparison of Computed and Experimental Forces & Moments. 


C L versus a 

a, 

degrees 

c 

^"Unumerical 

c 

^"L/experimental 

AC l = 

C - c 

^"L, numerical ^"L, experimental 

^^-'L / ^-'^experimental 

10.004 

1.806 

1.824 

-0.018 

-0.99% 

14.000 

2.078 

2.104 

-0.027 

-1.28% 

19.937 

2.447 

2.500 

-0.053 

-2.11% 

27.908 

2.790 

2.889 

-0.099 

-3.41% 

28.120 

2.796 

2.938 

-0.142 

-4.83% 

29.877 

2.845 

2.988 

-0.143 

-4.80% 

32.993 

2.867 

3.031 

-0.164 

-5.40% 


Note: C L experimental linearly interpolated from Table 3 (experimental data) for 


constant a. 


C D versus a 

a, 

degrees 

c 

^-D/numericaD 

c 

^-'D / experimental 

ac d = 

C - c 

V —D / numerical V -"D, experimental 

/ ^-'D / experimental 

10.004 

0.2613 

0.2783 

-0.0170 

-6.11% 

14.000 

0.3365 

0.3621 

-0.0256 

-7.08% 

19.937 

0.4629 

0.5030 

-0.0401 

-7.97% 

27.908 

0.6310 

0.7017 

-0.0706 

-10.06% 

28.120 

0.6346 

0.7058 

-0.0712 

-10.09% 

29.877 

0.6693 

0.7401 

-0.0708 

-9.57% 

32.993 

0.7238 

0.7917 

-0.0679 

-8.58% 


Note: C 


D, experimental 


linearly interpola 


ted from Table 3 (experimental data) for 


constant a. 
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Table 7. Concluded. 



Note: C D experimental linearly interpolated from Table 3 (experimental data) for 


constant C L . 



Note: C m experimental linearly interpolated from Table 3 (experimental data) for 


constant a. 
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Figure 1. Features of multi-element high-lift flows. (ML = 0.2, a ~ 34°) 
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Note: The origin is at the 
intersection of the leading edge 



Figure 2. Cutaway view of the computational domain and grid. 
(Selected mesh planes deleted for clarity.) 
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(a) Full wing-body (top view) 


(b) Wing-body juncture (bottom view) 


Figure 3. Experimental model installed in the NASA Langley 14x22 Foot 
Wind Tunnel. 
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Slat and flap in stowed (cruise) positons 

( Note: This configuration not investigated, 
for reporting wing dimensions only) 



Slat and flap in deployed (high-lift landing) positions 


Configuration Dimensions 

(Quoted with the slat and flap in the stowed (cruise) positions) 

Wing semi-span: b/2 

85.1 inches 

Mean aerodynamic chord: c 

39.6 inches 

Semi-span wing area: S 

3172.03 inches 2 

Aspect ratio: AR = b/ c 

4.56 

Leading edge sweep 

33.0° 

Quarter chord sweep 

30.0° 

Wing taper ratio 

0.4 

Slat chord 

5.0 inches 

Flap chord ratio: c J c local 

0.3 

Wing-root thickness ratio: (t/c) at root 

0.121 

Wing- tip thickness ratio: (t/ c) at tip 

0.1 

Body pod length (x direction) 

118.72 inches 

Body pod height (z direction) 

27.961 inches 

Body pod thickness (y direction) 

5.95 inches 

Leading edge of body pod (x,y,z) 

-17.70, 0.0, 0.0 (all in inches) 

Moment center (x,y,z) 

34.342, 0.0, 0.0 (all in inches) 

Slat bracket locations (experimental 

0.13, 0.33, 0.47, 0.64, 0.77, 0.94 (all in 

model only): y/(b/2) 

inches) 

Flap bracket locations (experimental 
model only): y/(b/2) 

0.13, 0.37, 0.61, 0.80 (all in inches) 


Note: The zero reference point for the configuration is defined as the intersection 
of the leading edge of the wing in the cruise configuration with the y = 0 plane. 

Figure 4. Model dimensions with sketches of a typical airfoil in both the cruise 
(stowed) and landing (deployed) configurations. 
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Figure 5. Airfoils at selected span stations with the slat and flap 

deployed (i.e., with the wing in the configuration for the 
present investigation). 
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Figure 6. Surface grid mesh. (Every other line deleted for clarity.) 
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(a) Complete wing. 



Figure 7. Typical grid mesh in the vicinity of the wing. (Vertical cut) 
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Method 

Grid 

Turb. 

Precision 



size 

model 


o 

TLNS3D 

12million 

SA 

Single 

□ 

TLNS3D 

12million 

SA 

Double 

A 

TLNS3D 

12million 

SST 

Single 

o 

TLNS3D 

1.5million 

SA 

Single 

> 

TLNS3D 

06million 

SA 

Double 




Figure 9. Force and moment plots. (For the experimental data, the balance accuracies in 
terms of aerodynamic coefficients were approximately: normal direction, + / - 
0.0024; axial direction, +/- 0.0017; and pitching moment, +/- 0.00095.) 
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Figure 9. Concluded. 
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Figure 10. Pressure coefficient distributions on the wing. (CFD 
calculations were made with the Spalart-Almaras 
turbulence model.) 
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(b) a = 14.000° 


Figure 10. Continued. 
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X 


19.937° 


Figure 10. Continued. 
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(d) a = 28.120° 


Figure 10. Continued. 
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(e) a = 29.877° 
Figure 10. Continued. 
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(f) a = 32.993° 
Figure 10. Continued. 
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(g) a = 33.980° 


Figure 10. Continued. 



(h) a = 34.858° 
Figure 10. Concluded. 
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Figure 11. Sketch depicting simulated oil flows and streamlines 
constrained to a vertical plane with the model at a 
very high angle of attack. (M^ = 0.2, a ~ 35°) 
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(a) Surface pressure-coefficient contours and constrained 
streamlines immediately above the surface. 


Figure 12. Model surface and flow-field parameters at a = 10.004°. 
(Spalart-Allmaras turbulence model) 
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Leading edge region 


Trailing edge region 


(b) Mach number contours and constrained streamline pattern for a stream- 
wise cut at n = 0.70. 


Figure 12. Continued. 





Leading edge region Trailing-edge flap 


(c) Velocity vectors for a stream-wise cut at ri = 0.70. 


Figure 12. Continued. 
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(d) Computed velocity profiles near the model surface at ri = 0.70. 


Figure 12. Concluded. 
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(a) Surface pressure-coefficient contours and constrained 
streamlines immediately above the surface. 


Figure 13. Model surface and flow-field parameters at a = 14.000°. 
(Spalart-Allmaras turbulence model) 
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M 


local 



Full wing 



Leading edge region 



Trailing edge region 


(b) Mach number contours and constrained streamline pattern for a stream- 
wise cut at r| = 0.70. 


Figure 13. Continued. 
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(c) Velocity vectors for a stream-wise cut at r| = 0.70. 


Figure 13. Continued. 
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(d) Computed velocity profiles near the model surface at r] = 0.70. 


Figure 13. Concluded. 
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(a) Surface pressure-coefficient contours and constrained 
streamlines immediately above the surface. 


Figure 14. Model surface and flow-field parameters at a = 19.937°. 
(Spalart-Allmaras turbulence model) 
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M local 



Full wing 



Leading edge region Trailing edge region 


(b) Mach number contours and constrained streamline pattern for a stream-wise 
cut at r| = 0.70. 


Figure 14. Continued. 
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Main wing element 



(c) Velocity vectors for a stream-wise cut at r\ = 0.70. 


Figure 14. Continued. 
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Phantom rakes (not physical) 


dT 



x = 0.50 
x - 0.85 
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v/V„ 



w/V. 


(d) Computed velocity profiles near the model surface at r] = 0.70. 


Figure 14. Concluded. 
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Top 




(a) Surface pressure-coefficient contours and constrained 
streamlines immediately above the surface. 


Figure 15. Model surface and flow-field parameters at a = 28.120°. 
(Spalart-Allmaras turbulence model) 


82 



M local 





Leading edge region 


Trailing edge region 


(b) Mach number contours and constrained streamline pattern for a stream- 
wise cut at r| = 0.70. 


Figure 15. Continued. 
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Main wing element 



(c) Velocity vectors for a stream-wise cut at r\ = 0.70. 


Figure 15. Continued. 


84 


Phantom rakes (not physical) 





x = 0.50 
x = 0.85 






v/V 


w/V. 


(d) Computed velocity profiles near the model surface at r\ = 0.70. 


Figure 15. Concluded. 
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(a) Surface pressure-coefficient contours and constrained 
streamlines immediately above the surface. 


Figure 16. Model surface and flow-field parameters at a = 29.877°. 
(Spalart-Allmaras turbulence model) 
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Leading edge region 


Trailing edge region 


(b) Mach number contours and constrained streamline pattern for 
stream-wise cut at ri = 0.70. 


Figure 16. Continued. 




Leading edge region 


Flap 


(c) Velocity vectors for a stream-wise cut at r| = 0.70. 


Figure 16. Continued. 
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Phantom rakes (not physical) 
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(d) Computed velocity profiles near the model surface at r| = 0.70. 


Figure 16. Concluded. 
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(a) Surface pressure-coefficient contours and constrained 
streamlines immediately above the surface. 


Figure 17. Model surface and flow-field parameters at a = 32.993°. 
(Spalart-Allmaras turbulence model) 
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Full wing 




Leading edge region 


Trailing edge region 


(b) Mach number contours and constrained streamline pattern at ri = 0. 


Figure 17. Continued. 
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Main wing element 



Leading edge region 


Flap 


(c) Velocity vectors for a stream-wise cut at r\ = 0.70. 

Figure 17. Continued. 
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Phantom rakes (not physical) 




x = 0.50 
x = 0.85 






v/V 


w/V 


(d) Computed velocity profiles near the model surface atr\ = 0.70. 


Figure 17. Concluded. 
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Top 




(a) Surface pressure-coefficient contours and constrained 
streamlines immediately above the surface. 


Figure 18. Model surface and flow-field parameters at a = 33.980°. 
(Spalart-Allmaras turbulence model) 
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(b) Mach number contours and constrained streamline pattern for a stream- 
wise cut at r| = 0.70. 


Figure 18. Continued. 
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Leading edge region 


Flap 


(c) Velocity vectors for a stream-wise cut at ri = 0.70. 


Figure 18. Continued 
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Phantom rakes (not physical) 
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v/V 


w/V 


(d) Computed velocity profiles near the model surface at r| = 0.70. 

Figure 18. Concluded. 
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(a) Surface pressure-coefficient contours and constrained 
streamlines immediately above the surface. 

Figure 19. Model surface and flow-field parameters at a = 34.858°. 
(Spalart-Allmaras turbulence model) 


98 



Full wing 




Trailing edge region 


(b) Mach number contours and constrained streamline pattern fora 
stream-wise cut at r\ = 0.70. 


Figure 19. Continued. 
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Main wing element 




Leading edge region 


Flap 


(c) Velocity vectors for a stream-wise cut at r\ = 0.70. 


Figure 19. Continued. 
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Phantom rakes (not physical) 




x = 0.50 


x = 0.85 
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(d) Computed velocity profiles near the model surface atri = 0.70. 


Figure 19. Concluded. 
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(a) a = 10.004° and 19.937° 


Figure 20. Effect of angle of attack on the trailing-edge flap. 
(Spalart-Allmaras turbulence model) 
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(b) a = 29.877° and 34.858° 


Figure 20. Concluded. 


29.877° 


34.858° 


103 



a = 10.004° 



a = 19.937° 


(a) a = 10.004° and 19.937° 


Figure 21. Effect of angle of attack at the juncture of the leading 
edge of the wing and the fuselage. (Spalart-Allmaras 
turbulence model) 
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a = 29.877° 



a = 34.858° 


(b) a = 29.877° and 34.858° 


Figure 21. Concluded. 
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(a) a = 10.004° and 29.877° 


Figure 22. Effect of angle of attack on the wing tip vortices. 
(Spalart-Allmaras turbulence model) 
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(b) a = 32.993° and 34.858° 


Figure 22. Concluded. 
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Figure 23. Effect of turbulence model on the wing 
pressure coefficients. ( a = 29.877°) 
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(a) Surface pressure-coefficient contours and constrained 
streamlines immediately above the surface. 


Figure 24. Effect of the turbulence model on the flow field in the vicinity 
of the flap, (a = 29.877°) 
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(b) Mach number contours and constrained streamlines for a stream- 
wise cut at r| = 0.70. 


Figure 24. Continued. 
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Spalart-Allmaras turbulence model 



Menter’s SST turbulence model 


(c) Velocity vectors for a stream-wise cut at r| = 0.70 
Figure 24. Concluded. 
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(a) r| = 0.70, x = 0.50 


Figure 25. Effect of the turbulence model on the velocity profiles 
near the wing surface, (a = 29.877°) 
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Figure 25. Continued. 


0.5 1.0 


113 


Phantom rake 
(not physical) 


Spalart - Allmaras trubulence model 
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Figure 25. Continued. 
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Figure 25. Concluded. 
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Spalart-Allmaras turbulence model 


Figure 26. Effect of the turbulence model on the wing 
vortices, (a = 29.877°) 
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